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ABSTRACT 


A formulation is proposed for the estimation of aerodynamic 
coefficients (parameters) of external stores through the 
technique of parameter estimation. The release of stores 
from the aircraft acts as an input and the resulting aircraft 
response is analysed through the Gauss-Newton method to 
estimate 'the store parameters. An example with simulated data 
has been used to show how the accuracy of estimation is 
affected by varying noise levels in the measured response, 
initial values of parameters used and the location of stores 
on the parent aircraft. Two possible approaches, called 
Hethod 1 and Method 2 are applied and compared : Method 1 
obtains the aircraft as well as store parameters from a 
single aircraft response due to store release; Method 2 'uses 
a two step approach - first aircraft parameters are obtained 
by usual input in the form of anelevator deflection and 
then only the store parameters are determined from the 
response due. to store release wherein the aircraft parameters 
are kept fixed at a priori values obtained in step 1. Method 1 
is shown to work well when measurement noise is either absent 
or its intensity is less while Method 2 is found to yield 
satisfactory results even in presence of high intensity 
noise in the measured data. 
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INTRODUCTION 
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Combat airplanes are normally equipped with a vast 
variety of external stores like bombs, fuel tanks, missiles 
and other ordnances to accomplish its operational missions. 
These stores are normally carried externally on wing and/or 
fuselage mounted pylons of carrier aircraft. No matter how 
potentially useful the store, if it damages or destroys the 
carrier or itself at release, it is of no practical value. 
Therefore, before any new store integration on an aircraft 
is cleared, it calls for detailed studies of handling 
qualities of the aircraft during carriage phase and also the 

separation characteristics of the store following its release. 

. 1 

Of many such studies reported in the; literature, covert 

and Schindel 2 have described the conditions for safe seperation 

of external stores. 

The reliability of predicted seperation characteristics 
depends on an accurate estimation of aerodynamic coefficients 
(parameters) of the store in the presence of interference 
flow field of the carrier aircraft and other external stores 
present in the vicinity^* Information on estimation methods 
for aerodynamic coefficients of stores of complex shapes and 
in presence of complicated interference flow field effects 
is scarcely available in the open literature# Theoretical 
methods/ currently available are restricted to steady spate 
conditions of the flow field around the model# Th^ prevailing 
wind tunnel methods such as captive model testing 


and 
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drop model technique have been used to generate aerodynamic 

data for the preflight simulation studies. Wind tunnel drop 

model, technique seems to be closer to the real situation 

than the captive model testing. However, both these methods 

are incapable of simulating the pronounced dynamic conditions 

of the stores encountered during release from the aircraft. 

Notwithstanding these limitations, a detailed survey of the 

methods and wind tunnel results on the aerodynamic loads 

associated with the external carriage of pylon-mounted 

stores adjacent to wing - fuselage combination has been 

reported by Marsden and Haines , These methods for both the 

subsonic and supersonic speeds are reviewed in Ref. 7, giving 

indication of their accuracy, range of validity and extent 

to which they have been verified experimentally. Further, 

variation of loads with store position, wing geometry, Mach 

0 

number is also discussed . Maddox et. al have compared 
flight results for captive loads with the corresponding 
results from several wind tunnel tests as well as with the 
most agreeable mathematical model when conditions were 
matched as closely as possible. In the above study, a store, 
similar in shape to Mk 83 bomb was mounted on a completely 
instrumented F-4 aircraft. The flight conditions spanned 
Mach 0.6-0. 9 in both maneuvering and steady flight. The data 
showed good correlation between flight test and wind tunnel 
results for moderate subsonic Mach numbers but pronounced 
divergence as the Mach number was increased. 
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It is, thus, clear that a need exists to obtain aerody- 
namic coefficients of the captive stores in real flight 
conditions. To that purpose, wa have proposed a method for 
estimating aerodynamic parameters of captive stores by 
employing parameter estimation technique of extracting 
parameters from input-output data of an aircraft. An excellent 
review of various methods used for parameter estimation from 
flight data is given by Maine and Iliff^. A brief outline 
and relative merits of various parameter estimation methods 
is given below. 

Aircraft parameter estimation is the process of deter- 
mining stability and control derivatives (Parameters) from 
the measured response of the aircraft for a known input. With 
the availability of high speed digital computers, parameter 
estimation has seen extensive practical applications. Following 
approaches are gene-'rally adopted for estimating stability 
and control derivatives : 

. Theoretical methods 

. Wind tunnel testing 

. Flight testing 

For tile purpose of preliminary design of an airplane/ 
store , theoretical methods arc useful, inspite of their 
limited accuracy. The wind tunnel testing improves the 
accuracy of estimation but it is a time consuming and cost 
prohibitive process of determining parameters. Accuro.be 
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simulation of control surfaces, power effects and flight 
conditions is difficult. Wind tunnel models used for most 
of the testing are often slightly different from the actual 
flight vehicle because of last minute configuration changes. 
Reynolds number differences and presence of support system are 
major reasons for discrepancies between flight and wind 
tunnel results. For these reasons, it is always desired 
that the wind tunnel estimates be corroborated with the 
estimates of actual flight testing. 

Various estimators available to extract parameters from 
flight data can be put under the following sub heads : 

. Equation error methods 

. Out put error methods 

. Advanced (Statistical) methods 

The-: equation error methods are based on the principle 
of least squares. It minimises square of the error in 
satisfying the equations of motion with respect to unknown 
parameters. Each equation is solved independently. The main 
appeal of the method is its computational simplicity and 
easy application to any linear or nonlinear model. The 
disadvantage is that the method can not be directly applied 
if all states are not measured accurately and produces 
poor results if the measurements are noisy. Considerable 
effort is, therefore, required for data reconstruction and 
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smoothing in order to obtain accurate results. These data 
preprocessing tasks are sometimes more complicated than the 
parameter estimation using equation error. However, these 
methods can be very effectively used as start up methods 
for more advanced estimators. 

The output error methods minimize the error between 
the measured and the model response produced for an identical 
input. It is assumed that the measured response is corrupted 
only by noise and that there are no modelling errors (process 
noise). Methods of this kind are capable of processing the 
measurement noise while assuming the model to be exact 
representation of the given system. These are probably the 
most widely used estimators for determining aircraft para- 
meters. A comprehensive survey of these methods is reported 

<| A 1112 

by Maine & Iliff . Newton-Raphscn method * with its 

9 13 

variants, gradient methods and Analog matching method 

fall in this category. A brief description of few of these 

methods is given in succeeding paragraphs. 

11 

Newton-Raphson optimization algorithm is the basis 
for most of the second order methods (methods that use 
second derivatives of the cost function) . It uses a two 
term expansion of the Taylor series. Since the evaluation 
of second gradient matrix is complex, the minimization of 
cost function is, generally very slow-. Moreover, if the . 
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second gradient is not positive definite, then the approxi- 
mating function does not have a unique minimum and the 
algorithm is likely to behave poorly. The performance of 
the method in close neighbourhood of a local minimum is, 
still, excellent. If the initial estimates are far from the 
minimum, the algorithm often converges erratically or even 
diverges. Further, there is necessity of inverting the 
second gradient matrix. The crucial issue concerning the 
inversion is that the matrix could be singular or ill 
conditioned. 

Modified Nowton-Raphson methods 11 are used where explicit 
evaluation of the second gradient of the cost function is 
complicated or costly but the performance of the Newton- 
Raphson algorithm is dasired. These methods approximate 
the second gradient in terms of the first gradient; The 
approximation generally improves the speed of convergence. 

Advance methods are capable of determining parameters 
in the presence of measurement and/or process noise. These 
methods are based on the probabilistic concepts. Maximum 
a posteriori probability (MAP) estimator and Maximum 
likelihood (ML) estimator 9 ' 14 belong to this category. 

The latter estimator is one of the most favoured estimator 
because it yields parameter estimates that arc asymptotically 
unbiased, consistent and efficient. Whan process noise is 
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absent and the covarience of the measurement noise is known, 

ML methods reduce to output error methods. In absence of 
measurement noise, ML methods reduce to Equation error 
methods. The algorithms accounting for both the measurement 
as well as process noise require more computer time and 
core. Moreover, convergence problems and other practical 
difficulties are often encountered. 

< C -1 fZ 

Recently, Raisinghani and Adak ' have proposed a 
simplified output error method, called Gauss-Newton (GN) 
method for aircraft parameter estimation. In this method, the 
model response is expanded in a Taylor series about the 
current trial value of the parameters and retains only the 
linear terms. This linearized model is substituted in the 
least square objective function and solved in the least 
square sense. It is shown that the method works satisfactorily 
even for the noisy data provided 'the initial trial values 
of the parameters are not too far off(upto + 50 % off true 
values). Since the computational effort involved is consi- 
derably less than the more advance methods like maximum 
likelihood m ethod, while the accuracy obtained is fairly 
good, it was decided to use this method as 'the parameter 
estimator for the present work. 

The problem studied here can be divided into two 
distinct phases : (i) to postulate a mathematical model. 
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(ii) to use a suitable output-error or any other estimation 
technique to extract aircraft and/or store parameters. The 
problem has been formulated by treating release of stores 
as the step input which excites the airplane dynamics to 
provide the output in the form of aircraft motion variables. 

The magnitude and direction of perturbations will be a 
function of the airplane's dynamic characteristics, type of 
stores, location of stores and the release/jettisoning 
mechanism. The equations of motion of the aircraft for such 
an input have been derived and include store parameters, 
in addition to aircraft parameters, as unknowns. Only 
longitudinal parameters of the aircraft and store ware 
estimated. It is shown that an adequate excitation of aircraft 
dynamics by store dropping will enable simultaneous estimation 
of store and aircraft parameters. However, for locations 
and/or type of stores leading to poor excitation of aircraft 
dynamics and high noise levels, two step approach is 
recommended : (i) estimation of aircraft (without any stores 
attached) parameters through control (elevator) input 
(ii) store parameters estimation from aircraft response due 
to store dropping wherein all the aircraft parameters are 
fixed at values estimated from the first stop. This approach 
was also used to increase the accuracy of estimation when - 
eve r the noise level in the measured data was high. 
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The problem of aircraft response following store release 

17 18 

was formulated by Rao and Raisinghani and Rao . In this 
formulation, drag coefficient, lift coefficient and pitching 
moment coefficient of the store appear alone or in combi- 
nation as equivalent of control derivatives in the equations 
of motion of the aircraft. In the present work, a more 
direct approach has been used to arrive at essentially 
the same equations of motion except for correcting a few 
errors and inconsistencies found between Refs. 17 and 18. 

It may be emphasized that the methodology proposed for 
store parameter estimation can be used in conjunction with 
any one of the methods available in the literature for aircraft 
parameter estimation. . 

Duo to nonavailability of real flight data, the 
proposed methodology for store parameter estimation has been 
validated on simulated data only. An aircraft, resembling 
FIAT-G91 and missile similar to those reported in Ref. 7 
were used for the case study. Pour different locations of 
stores were used to simulate flight data following stores 
release. Simulated flight data were corrupted with pseudo 
noise of various intensities. Store parameters were extracted 
either along with the aircraft parameters or alone while 
fixing the aircraft parameters at values estimated separately. 
Details of the effect of such one step or two step estimation 
process arc reportad for various locationsand for different 
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intensities of noise. For most of the cases, it is shown 

that there is a close agreement between the estimated and 

true values of the parameters and also the Cramer-Rao bounds 

deviation 

(°"Cr) are ^ ow * M; - an values and sample standard /( <5 ~ s ) of 
parameter estimates arc also estimated and compared with 
Cramer-Rao bounds ( < 3 - for various noise levels. The 
agreement between r and is found to be good. Finally, 

it is shown that the two step- procedure, is able to estimate 


the store parameters quite accurately even when the . initial 
values of the store parameters were set arbitrarily as f-ar 
of as + 500% off the true values and also in presence of 
high intensity noise. 

The formulation of the equations of motion of aircraft 
following store release is given in Chapter 2. Details of 
the Gauss-Newton method used for parameter estimation are 
also given in same Chapter. Chapter 3 contains the details 
of the application of the method to a test case. The results 
and discussion are given in Chapter 4. Main conclusions of 
the present study are listed in Chapter 5. 
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CHAPTER - II 
FORMULATION 
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The estimation of stability and control parameters from 
flight data postulates a mathematical model of the aircraft 
under test. An important question relates to the complexity 
of the assumed model. Although, increasing the complexity 
of the model (i.e., increasing the number of unknown para- 
meters) may lead to better description of the aircraft motion, 
it may result in too many parameters being sought from a 
limited amount of data and thereby lead to reduced accuracy 
of estimates. It would be, therefore, our endeavour to 
postulate an appropriate model for the problem to be studied 
here . 

For the present study, the perturbed motion of the aircraft 
following the stores release can be expected to be small. 
Further, we will be interested in analysing the aircraft 
response of a short duration immediately following the 
stores dropping. The airplane with external stores (called 
‘Loaded aircraft' here after) is assumed to be rigid and in 
a steady state, rectilinear, wing level flight. Thus, we . 
assume that the perturbed motion of aircraft due to stores 
dropping can be represented by the usual decoupled set of 
equations of motion for longitudinal and lateral-directional 
motions. Since we would consider the case of symmetric load 
dropping only, the perturbed motion will be considered only 
in the plane of symmetry. In calm atmosphere, the longitudinal 
perturbed equation of motion in the stability axes system 
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can be written as follows 

■ x(t) = A x(t) + B q (t) 
where 


x(t) 

= fu 

oC ^ . 6 

- T 

f 

• 

x(t) 

= fu 

ci q o 

] T 


>3(t) 

= [«e 

i 

i 



i 

1 

fx 

u 

X, 

O' 

0 

-gcose^ 


Vb 

! . 

Zj/Vx 

1+s i 

-gs ine^/U-L 

. : 

■ -■ . i 

I 

M u 

My 

M q 

0 

' ' ; 

. | 

1 0 

0 

l 

0 

= 

r 

1 X r 

L 

A ^1 

~ e 

M . 

£e 

- T 

0 ; 1 

' 5 


( 2 . 1 ) 


and 

B 


Matrix, x contains th© state variables, ">> (t) is the 
control input, matrices A and B contain the dimensional 
stability and control derivatives as defined in the nomen- 
clature. 

The algebraic/transcendental equation of measurement can 
be written as 


y (t) = h 


x(t), 'n(t) 


, | 


( 2 . 2 ) 
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where y(t) is the vector of measurements, h T.~l defines 

l- _J 

■the relationship between the measurements, the state and 
the control, and *b(t) is a vector of additive measurement 
noise. Although all or some of the output variables may 
be different from the state variables, for the present 
study using simulated data only, we will assume for simplicity 
that either all the state variables are measured directly 
or have been obtained through a state estimation algorithm. 
Thus , the measurement equation is written as 


y^ (t) = x^(t) + uk (t) (2.3) 

where y^(t) and x^(t) are the components of measurement 
and state vectors and ui^(t) are the measurement noise 
components . 


It must be pointed out that for all the output error 

methods, including the present one, the modelling error 

19 20 

must also be included as measurement noise * . For real 

flight data, this may tend to make the noise quite coloured. 
However, for simplicity, we have assumed that the measurement 
noise components are independent, stationary Gaussian noise 
sequences . Again, in real data, there may be biases in the 
flight data, which would have to be included in Eq. (2.3) 

and estimated through the estimation algorithm. For simplicity 

; ' - ’ ■ 

again, we assume that our simulated data have no biases. 
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The above set of equations of motion were used in 
Refs. 15,18/21 to estimate aircraft parameters from- flight 
data. A known control (elevator) input was applied to obtain 
the perturbed response of the aircraft and Gauss-Newton 
method was used to extract parameters from- it. For the stores 
release case, the controls are assumed to be locked while the 
stores are dropped. Thus, the input to the aircraft is in 
the form of a step input due to dropping of the stores. 

The equations of motion corresponding to such an input were 
obtained by Raisinghani and Rao . However, we shall obtain 
the same set of equations in a more straight forward manner 
by considering the input forces and moments acting on the 
airplane following the stores release. The following changes 
take place when stroes are released symmetrically from an 
aircraft. 

(1) Gross weight reduction equal to the weight of the stores 
dropped. 

(2) Centre of gravity (C.G.) of aircraft shifts, both 
vertically and longitudinally. 

(3) Moments of inertia change. 

(4) Contribution of forces and moments to the loaded aircraft 
due to stores ceases to act following stores release. 

The magnitude of those forces and moments depends upon 
the type of stores, location of stores, aircraft 
geometry, shift in C.G. and flight conditions etc. 
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(5) The area of the wing from where the stores are dropped 
becomes aerodynamically cleaner and, thus, changes forces 
and moments contribution of the wing. 

(6) Till stores remain in the immediate vicinity of the wing 
after their seperation, interference effects cause 
variation in forces and moments acting on the wing. 

Items 5 and 6 in the above list are very complex and 
difficult to estimate theoretically with any degree of 
confidence. It is hoped, however, that change in forces and 
moments due to these effects will be much smaller as compared 
to other effects listed above. Therefore, these have not been 
included in the formulation. 

Following release of stores, there will be a shift in the 
C.G. location. The magnitude and direction of shift in C.G. 
will be a function of the type of stores and their location 
on the aircraft prior to release. For symmetric release of 
stores, C.G. of aircraft is assumed to shift only in the 
plane of symmetry. Since stores are mostly so located that 
the longitudinal C.G. locations of stores and of the aircraft 
are almost coincident^ the change in the longitudinal location 
of C.G. following storbs dropping is negligible and, therefore, 
not considered in the present formulation. However, the 
vertical shift in C*G. location following stores release is 
considered. It may be mentioned that the longitudinal shift 
in C.G. can be easily included in the formulation, if 




desired for a particular configuration. 



Let us consider an even number of identical stores 
located symmetrically under the wing of an aircraft. The 
steady state flight of such a configuration will be consi- 
dered first. The equilibrium condition implifcs that the net 
forces along X and Z axes, and the pitching moment about 
Y-axis is zero. These forces and moments are due to the thrust, 
weight and aerodynamic forces of the- aircraft and stores. 

Let the drag, lift and pitching moment acting at C.G. of the 

store be denoted by D , L and M respectively. The distances 

x s s s 

between the store and aircraft C.G. are shown in Pig. 1. The 
changeis in the drag, lift and pitching moment of the aircraft 
due to store dropping can be written as follows : 
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(2.4c) 


where N is the number of store (2 in our case), m is the 

S o 

mass of the store, X is the distance of store C.G. from 
the aircraft C.G. (positive behind), Z is the distance of 


store C.Gi below the aircraft C.G. (positive 'downward), 
d T is the vertical shift in the aircraft C.G. (positive 
gteWQv fS.£e Fig. 1). The last term on the right - 
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hand side of Eqi (2.4c) arises due to the vertical shift of 
aircraft C.G. following the stores release; its occurance 
is explained as follows. 


Prior to release of stores, the thrust and drag forces 

contributing to pitching moment were equal and opposite. 

After release of stores the net drag force is reduced by 

D N . Equivalently, we could say that the thrust force 
S s ~ 

exceeds drag force by an amount equal to D N * Due to the 

s s 

vertical shift (d„) in C.G. of aircraft, this excess thrust 
will contribute additional pitching moment about C.G. given 
by the last term of Eq. (2.4c). 

The . above changes in forces and moments following stores 
release are treated as step inputs to excite the aircraft 
dynamics. Thus, we can treat these inputs equivalent to a 
step control inputs whose magnitudes are determined as 
follows 
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% S C C M. 


M s. 


""e 




yy = 


yy 

J 2 -a. ^'e 


(2.5c) 


As defined above, the control derivatives X r and Z f 

0 0 

represent for a unit step change in £ , the change in forces 
per unit mass of aircraft along X-axis and Z -axis respecti- 
vely. Similarly, M, represents for a unit step change in > 

* 0 

the change in pitching moment per unit I . Using the 
changes in forces and pitching moment arising duo to the 
stores release, the equivalent control derivatives for the 
present problem can be defined as follows . 
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(2.6c) 


Thus 


, A g and will replace respectively X g , and 


M r in Eq. (2.1) with v | (t) replaced byl(unit step) . 

Alternatively, the above, expressions could be arrived at 

17 

by following the formulation similar to that of Rao 



The drag, lift and pitching moment of store can be 
expressed in terms of nondimensional coefficients as follows 


D S = C D S % A F 


L S C L„ cm S_ 
S H S S 


M S ~ C M q S S S C S 
o 


( 2 . /a) 

(2.7b) 

(2.7c) 


where q c is dynamic pressure seen on the store and is given 

O 

by 


q S = q l S CF 

Here q 1 is steady state dynamic pressure 
factor whose value is a function of type 
of store and the Mach number. 


and is 

of store. 


( 2 . 8 ) 

scale 

location 


a is frontal cross-sectional area of the store 
F 

S c is exposed area of two wings of store. 

O 

and C„ is mean chord of store wing. 

s 


Substituting Eq* (2.7) 

q s 

X, = Ay = 

Op " 


in Eq. (2.6) we get 
+ m s g Sin e^) 


(C f q fl S s - m s g cos e^) Ng 



( 2 . 9h ) 

(2.9b) 
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( “ C M = C S + C L V % S S N S +C D % A F W'V 

s s s 


YY- 


(2.9c) 


Once the above model of the aircraft is assumed to be 
known/ the system identification problem reduces to that 
of parameter estimation. As mentioned earlier/ the Gauss- 
Newton algorithm for parameter estimation proposed by 

•1 CZ A k" 

Raisinghani and Adak 1 * has been used for the present work. 
A brief outline of the method is presented hare for the sake 
of completion. 


Let the dimensional stability derivatives appearing in 

matrix A and control derivatives in matrix B be denoted by 

C, k=l,2, ...., m where m is the total number of unknown 

parameters of airepaft and store. State variables : u, , q 

and © are represented by x^(t), 1=1,2, ....,n where n is 

total number of state variables. Let the corresponding 

measured state variables be denoted by x . (t) . For the some 

mi 

input (elevator or store dropping) and using initial starting 

values of C k = C£ , the estimated response of aircraft is 

obtained by solving Eq. (2.1). Let this be denoted by x°(t) . 

The difference between the measured, x . (t) and estimated, 

mi 


O 

> 


,°(t) responses is attributed to the difference between the 


values of and C°„ The aim is to change C“ values such 

that x'i (t) — > x . (t) . Let the values of C? be changed by 
1 mi k 
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so that C k = C° The corresponding estimated 

.response for = C^. is denoted by x^ (t) . Following Gauss- 
Newton minimization procedure, the model is linearized by 
expanding x^(t) in a Taylor series about the current trial 
values of the parameters and retaining the linear terms only, 
we ge t 


(t) 


x9(t) + 


^x^(t) 


■>c 


r> 


. + 


*x ± (t) 

~bc r 

m * 


* C 1 + 


^x ± (tT 

TcT“ 


AC 2 + 


<aC 


m 


(2.10) 


Here superscript ‘o' means quantities evaluated at initial 
trial value. 


A least square objective function is written as follows 
n N— 1 p * 2 

minimize S = T V x . (t .) - x .(t .) j (2.11) 

t I mi j i j I 

i=l >0 <- 

where N represents the data points used for analysis of 
each of n state variables. 


Now the linearized model, Eq. (2.10) is substituted into 
the objective function, Sq. (2.11) and the equations are 
formed by setting the partial derivatives of S with respect 
to each of the unknown parameter equal to zero. 






0 , k — 1,2,...., 
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The resulting equations for the longitudinal case are as 
follows 


(P T P) C = P T ( a X) 


( 2 . 12 ) 


where superscript 1 T 1 designates the transpose of the 
matrix and 
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AC = 


AC 


1 


AC. 


T 


aC 


m 


I 

i 


u m (t o ) " u ° ( ’ t o ) 


u m ■ 

* - (t o } 


A X 


%n (t o? " q ° (t o ) 

q n/ t N-P “ q ° ^N-p 

s m (t o ) - e ° (t o ) 


“mN-l 5 - e ' (t N-l > 


The partial derivatives in matrix, P are function of 
time and parameters, C£, , To solve the above set of Eqs.(2.12) 
for aC^'s, we need these partial derivatives. The following 
procedure is used to obtain the desired partial derivatives as 
a function of time. 
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The governing differential Eq. (2.1) is differentiated 
with respect to each of the parameters being estimated. . 

First terms on the left-hand side of these equations will 
be of the type ) whose order of differentiation 

is interchanged. On the right hand side/ the terms representing 
variation of input with parameters/ C, (such as etc) are 

K c A 

U 

set to zero as the control input is independent of the changes 
in parameter. Thus, we obtain the following set of 12 equations 
of the form similar to Eq. (2.1), i.e. 
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with r = u, y , q and used to obtain the required 12 sets 

of equations which can be solved to obtain the 48 partial 
derivatives and these, are used to form the matrix P in 
Eq. (2.12). 


Now the set of linear algebraic Eq. (2.12) can be 
by any appropriate technique for AC, 1 s. These are use 
update the current estimates for by writting 


solved 
d to 


C k I+1 = +aC ] c ' , I = 1/2/ ... . / MI 

where NI being the number of iterations. 

Using the so obtained updated values of parameters/ the 
aircraft response is calculated and matched with measured 
response . This process is continued until the matching 
between the estimated and measured flight response is within 
the specified margin or the number of iterations exceed the 
specified limit/ whichever occurs earlier. 


Since the partial derivatives arc evaluated with respect 

to non-optimal parameters/ C^, at each iteration, fhe parameter 

improvement, does not immediately lead to the optimal 

values of C ‘ s . However, if the process is convergent, the 
K 

&C improvement will tend to zero such the t the es time tec 
k 

airc raft response will tend towards the measured response 
and the ck 1 values will approach the optimal values as the 

K 

number of iterations increases. 



CHAPTER 


III 


APPLICATION OF THE METHOD AND CASE STUDY 
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3.1 TEST CASE 

The method formulated in previous chapter was applied 
on a test case. An aircraft resembling FIAT G-91 was 
selected. Its inertia and geometric characteristics and 
stability and control derivatives 13 are listed in Table 1. 

The store was selected from Ref. 7 so that its aerodynamic 
coefficients could be estimated using the graphs given in 
that reference. The geometric details of the store are also 
given in Fig. 1, 

Since/ no real flight data could be obtained for the 
aircraft response following stores release, the method was 
applied on the simulated data. For this purpose, a computer 
program in Fortran IV was developed. A flow chart of the 
computer program used for parameter estimation is shown in 
Fig. 2. The main features of the program are as follows : 

(i) For the true values of the parameters and known input 
form, it generates simulated measured response of the 
aircraft. The input can be given either in the form of an 
arbitrary elevator input (Table 2 and Fig. 3) or in the 

form of stores release. This kind of response will be referred 
to as no-noise response. 

(ii) Measurement noise of various intensities could be 
added to the no-noise response to generate noisy responses. 
Pseudo random numbers were computer generated so as to have 



30 . 


a normal distribution with zero mean and assigned standard 
deviation. The intensity of noise was varied to correspond 
approximately to 1%, 2 %, 5 % and 10% of the maximum magnitude 
of the. corresponding motion variable' . A typical set of 
standard deviation values (cr N ) of noise used for all the 
four store locations corresponding to 5% noise are shown in 
Table 3 . 

(iii) Using the initial values of the parameters assigned for 
the model, it generates model response by solving Eq. (2.1) 
for the same input as was used to generate the measured 
response. The input could be an elevator deflection or a 
step input due to store release. 

F* 

(iv) It calculates the partial derivatives (sensitivity 
coefficients) as required to form matrix, P of Eq. (2.12) 
by solving Eq. (2.13). 

(v) The set of Eqs. (2.12) .is arranged and solved to obtain 
the increaments, nC^s, for the parameters being estimated 
and the values of the parameters are updated. It has the 
option of utilizing either no-noise response or noisy 
response obtained in steps (i) and (ii) above. 

(vi) Using incremented parameters as initial values, steps 
(iii) to (v) are repeated each time, till the model response 
and measured response match within the specified accuracy 
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or the change in successive values of the parameters is 
less than the specified or the number of iteration exceeds 
the assigned limit. 

(vii) There is an option for using the above basic schema 
in the following two ways . 

(a) To estimate all the aircraft as well as store 
parameters from the measured response due to stores 
dropping in one step. This will be referred to as Method 1. 

(b) The parameters of the aircraft can be first estimated 
by analysing aircraft response to known elevator input. Next, 
these values of aircraft parameters are kept fixed at 
estimated values and the estimation algorithm used again to 
estimate only the store parameters from the aircraft response ' 
due to stores dropping. This we shall refer to as method 2. 

For solving Eq. (2.1) to obtain measured and model 
response, a fourth order Runge-Kutta method was employed. 

A step size of 0.005 second was used and the data was 
scrambled back to a step size of 0.1 second for the purpose 
of GN method. This procedure was also followed for obtaining 
sensitivity coefficients at step size of 0.1 second. A 
signal length of 4.9 seconds (50 data points) was generated. 

The pseudo numbers for noise contamination were generated 
by a built-in computer subroutine, GGNOR. The noisy measured 



responses with 1% to 10% noise levels were analysed by using 
both Method 1 and Method 2, The noise samples for all the 
motion variables were independent of each other. Further, for 
studying the effect of noise on parameter estimation, 20 
samples of noise with same mean and standard deviation were 
generated and added to the same no-noise response to prepare 
twenty samples of measured noisy responses. These twenty 
responses yielded twenty different estimates of parameters 
which were utilized to obtain statistical properties of the 
estimates. Specifically, the sample mean and sample standard 
deviation of parameter estimates was obtained as follows 


5 = 


1 

N 


N 

i=l 


(3-1) 


1 _ 

N- 


N 

7Z (C.-C) 
i=l 1 


h 


(3.2) 


where 

C = Mean of parameter estimates 
C.= Parameter estimate of ith sample 

l 

N = Number of parameter estimates (20 in present study) 

= sample standard deviation 

The computer program also contains provision for obtaining 
Cramer-Rao (CR) bounds fop each of the estimated parameter. 

The CR bound 12 provides the minimum variance for parameter 
with which it can be estimated for a given set of data. 
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Thus, it provides the amount of confidence to be placed in 
the various parameter estimates. Since we are dealing with 
simulated data where the standard deviation of noise (°“N) 
is known, the CR bounds could be obtained in the following 
two -ways . 


(i) Using the magnitude of the difference between the 


measured and estimated response, x - x the CR bound can 

m . €3 » 


be estimated as follows. 
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Where D is a diagonal matrix of size N x N with diagonal 
elements (d . . ) given by 


d. . = 
xi 


•N 

1=1 


x (t . )— 
m x 


( t ± ) j 2 /n 


1, N 


and x = u/U^, ^ , q, © 


(ii) Using the knowno- N for all the measured variables, the 
lower CR bounds in terms of variance for each of the 
parameter estimates is given by 



P T B 



(3.4) 
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where x = u/ ,<*■, q and © 

1 

Matrix P is given by Eq. (2.12a) 

Thus, it was possible to compare c-g with CR bounds, 
a- and ^- CR obtained by using Eqs . (3.3) and (3.4). 


3.2 NUMERICAL EXAMPI+E 

To validate the computer code, measured response of the 
aircraft following stores dropping was obtained for the 
various store locations given in Table 5. The parameters for 
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for the store (Fig, 1 ) were estimated from Ref. 7 ''for. the purpose 

of present study. The true values of store parameters: C t 

S 

and C M so obtained (Table 5) were used to calculate 

s s 

A.., A and A' as required to solve Eq. (2.1). Four store 
^ , V2C1 

locations, named V1C1, VlC2/and V2C2 were chosen. Notation 

VI & V2 refer to store C.G. lying vertically below aircraft 

C.G. by a distance of 0.15 C and 0.4 C respectively. While 

Cl and C2 refer to store C.G. lying longitudinally away 

from aircraft C.G. by a distance of - 0.25C and 0.25 C 

respectively. Store C.G. lying below and behind aircraft 

C.G. is taken positive. The 'spanwise location for all the 
fixed 

cases was half semi span position of the aircraft wing. 

To illustrate, details of numerical calculations employed 
for obtaining A x , A and for one particular location V2C2 are 
given below. This location means (refer to Fig. 1) 

X c = 0.25 C = 0.5125 m 

Y g = 0.5 (b/2) =2.15 m 

= 0.4 C = 0.82 m 

s 

Store drag coefficient, C D is assumed to be constant 

s 

and equal to 0.2 based on the frontal cross-sec tional area, 

A . The dynamic pressure (<3g) seen on the store is calculated 
from q g = c^S^,, where S^ F is obtained from Fig. 5 of Ref. 7. 

N is the number of stores dropped. For the store under 
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study we have 


nig = 5 0 kg 


Cg = 0.24 m 


; Ap = TL- d 2 = 0.0314 m 2 


; Sg = 0.23 m 


X CGS = 0.25 Cg= 0.06m ; N = 2 


= 0.0 rad 


; m = 5000 kg (refer Table 1) 

O, 


? = 0.685 kg/m" 


; S CF = 0.915 


U = 137.5 m/sec 


q S q l S CF 


? I = 30400 kg-m" 
% ¥ 


0,5? U 2 S cp = 5925.0 N/m' 


= 


( % q s a f + m s g sin @ 1 } “S 


N, 


= 0.01488 m/sec 


2 


m 


The values of the store lift coefficient/ C T _ (based 

JuD 

on the exposed area (S_). of two wings of the store) and the 

O 

store pitching moment/ (based on the mean chord (Cg) of 

s 

the store) obtained from Fig. 11 of Ref. 7 are : 


C = 0.7 , C = -0.2 

b b 

Now, A is calculated from the following equation 
Zj 

(C L = q S S S * m S 9 COS V N S , 

A,. = - — § — 0.1854 m/sec 

4 ™ 
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CHAPTER - IV 
RESULTS AND DISCUSSIONS 
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4 . 1 GENERAL 


In this chapter, we shall discuss the results of the 
test case for estimation of store parameters. The response 
of test airplane following stores release is utilized to 
estimate desired parameters through GN method. The airplane 
and the stores used for the study, the governing equations 
of motion and GN method were described in earlier chapters. 

The proposed method has been used to extract longitudinal 
store parameters from no-noise response as well as from noisy 
response of aircraft following stores release. Both Method 1 
and Method 2 described in Chapter 3 have been extensively 
applied on simulated noisy response and relative merits of 
these methods are pointed out. Attention was also focussed on 
the effect of noise level present in the measured response 
and the effect of initial values of parameters on the accuracy 
of the estimated parameters. To this purpose, the results 


obtained for noisy responses will be discussed under the 

following subheads. All the results using Method 1 are 

subhead 

discussed under^ mmm of case 1 whereas the results using 

subheads 

Method 2 are further divided under/ mmm of case 2 to 

case 5. These cases differ from each other in the way the 
two steps of method 2 are used, i.e., the difference in the 
approach used for the estimation of aircraft parameters in 
step one and/or the estimation of store parameters in 
step 2. For ready reference, approach used for cases 2 to 5 


are defined below 
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. Case 2 


. Case 3 


. Case 4 


. Case 


One measured noisy response fo-ran arbitrary 
elevator input was used for aircraft parameters 
estimation. The aircraft parameters were kept 
fixed at these estimated values while estimating 
store parameters from aircraft response due to 
stores release. 

Step 1 was again similar to case 2, i.e-. , aircraft 
parameters were estimated from one noisy response 
for an arbitrary elevator input. The aircraft 
parameters were kept fixed at these estimated 
values while estimating the store parameters from 
20 different noisy responses due to stores release. 

The aircraft parameters were estimated from 20 
different samples of noisy responses for an arbitrary- 
elevator input. The mean values of the aircraft 
parameters so estimated were used as the initial 
valups and kept fixed while estimating the store 
parameters from one noisy response due to stores 
release. 

In this case, step one of Method 2 was replaced 
by choosing true values of aircraft parameters and 
keeping them fixed while estimating the store 
parameters from one noisy response due to stores 
release. 
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The Cramer— 'Rao bounds («r^) for the estimated parameters 
were also obtained for cases 1,2,4 and 5. The mean and the 
sample standard deviation ( <r ~g) of parameter estimates for 
store parameters were Obtained for case 3. A comparison of 
a _ s andg-^ is also presented. 

The initial values of the aircraft parameters were 
arbitraily chosen to be +50% off the true values while the 
store parameters were +500% off the true values. Many such 
combinations were tested to validate the computer code; 
ho wsver, results will be presented for one such set of initial 
values shown in Table 4 and 5 : both the true values and the 
initial values of store parameters at different locations of 
the stores are shown in Table 5 . For convenience, the response 
of the model with the estimated values of the parameters will 
be referred to as the ‘estimate^ response* while that with 
the initial values of the parameters will be referred to as 
the ' initial response 1 . 

4.2 RESULTS FOR NO-NOISE CASE 

Aircraft and store parameters were estimated from no-noise 
measured response for all the four store locations. In all 
cases, parameter estimates converged very close to their 
true values in about five iterations. Since the estimated 
values and true values of all the aircraft as well as store 
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parameters matched up to four decimal places, no comparison 
of these Values is shown, in a tabular form. However/ a 
comparison of the measured and estimated response for all the 
four measured variables : u, o< / q and e for one typical 
location V1C1 is shown in Pig. 4. Initial response is also 
shown on the same figure for comparison. It is seen that the 
matching between the measured response and the estimated 
response is almost perfect. Thus# we may conclude that the 
proposed method can extract all the aircraft as well as 
store parameters quite accurately through GN. method (Via 
Method 1) from' no-noise response of aircraft due to stores 
release. 

4.3 RESULTS FOR NOISY CASES 

To study the effect of measurement noise on the accuracy 
of parameter estimation, noisy responses contaminated with 
1%, 2%, 5% and 10% noise were analysed. For convenience, 
results from noisy responses are presented separately for 
Method 1 and Method 2. A comparison between results for 
various cases will also be presented, where appropriate. 

4.3.1 Results for Method 1 and Case 1 

' -n’rrn,- — -• ."U in "'Tr.'-r,' • 1 ~ 1 ' ■ — 

Case 1 was studied for all store locations and for noise 
levels of 1%, 2% and 5%. The estimated values of the parameters 
and their CR bounds, both ctnd given by Eqs. (3.3) c:.nct 
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(3.4) for location V2C2 are shown in Table 6a, b. Table 6a 

shows that for low noise levels (up to 2%), all the parameters 

are well estimated exc-opt for the so called weak derivative, 

2 S 21 

' . Notwithstanding poor estimation of z , the estimated 

and measured responses matched quite closely as shown in 

Fig. 5 for one typical location V2C2 and for noise level of 

2%. It may be noted from Table 6a that the true values lie 

within + o- CR of the estimated values. Further, ^CR are quite 

small for low noise levels and, as expected, increases 

proportionately for higher noise levels. However, the estimated 

values of the store parameters C , C T and C for 5% noise 

"S S 

do not compare well with the true values and have relatively 
large Better results were obtained by use of Method 2 

as discussed in subsequent subsections. 

Cramer-Rao bounds defined as cr- CR and f'v'cR ^° r noase 

levels were also compared as shown in Table 6b. As expected, 
for the simulated data, the values of and are quite 

close to each: other. It was further desired that both cr- CR 
o-'nd (F- be compared with sample standard, deviation of the 

parameter estimates. For this purpose, one typical response 
data was contaminated with 20 different samples of noise 
having the same mean and standard deviation (c~^) • Analysing 
these 20 samples, the mean values of the parameter estimates 
and cr were calculated using Eqs . (3.1) and (3.2) respectively. 

s 
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Values so obtained are given in Table 7 for noise levels 1%, 

2% and 5%. It is noted that the mean values of all the 
parameters, particularly of the store 'parameters show consi- 
derable improvement. Further, the sample standard deviation 
also are smaller. 

A comparison among c~ Q , -~~ CR and 5- Cp for the estimated 

parameters for different noise levels is shown in Table 8. 

.7“ s ~c 

It shows that the ratios — - and ■=; — lie in the range 

° CR ^CR 

of 0.75-1.3 which is considered to be a reasonably good 
20 

agreement for a statistical sample of 20 only. 

The above results of Case 1 show that the store parameters 

estimation is not so accurate for high noise levels. It was 

*$ 

conjectured that this was due to an attempt to estimate too 
many parameters from insufficient information contained in 
the measured response being analysed. This ■ motivated the 
introduction of Method 2 for which the results are presented 
next. 

4.3.2 Application of Method 2 

To improve the accuracy of estimated store parameters, 
a two step approach, called Method 2 and explained in Chapter 3, 
v/as adopted. Airplane without stores is excited by an arbitrary 
elevator input (Table 2 and Fig. 3) to obtain airplane response, 
which is used to obtain airplane parameters through GN method. 
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The aircraft parameters and CR bounds so obtained for no-noise 
case and for noise levels of 1% to 10% are presented in 
Tables 9a and 9b. These are the values of aircraft parameters 
that were used as a priori fixed values in Case 2 and Case 3 
when step 2 was employed to estimate store parameters. Also, 
twenty different noisy responses for an identical elevator 
input were generated by adding 20 different noise samples 
of same mean and standard deviation ( cr-^) . These were also 
analysed through GN method to obtain mean values for aircraft 
parameters. These values are presented in Table 10 and 
were used as a priori fixed values in Case 4. From Tables 
9 and 10, it may be observed that the .so called weak parameter. 


Z is again poorly estimated, as was also observed earlier for 

Method 1. Further, another weak parameter, X c shows relatively 

® e 

poor accuracy of estimation. However, it may be pointed out 

that the value of X. is not required for step 2 of Method 2. 

° e 

It may be pointed out that another weak derivative Z ^ was 
estimated through M by relating and h Y 

means of the known geometry of the aircraft as given below: 


M 


s 


YY, 


1. m_ 
t a 


The weak parameters X . and Z are so called, since 

be H 

they do not affect the response significantly. That is why, 
inspite of poor estimates for ^ and Z q , the matching 
between the measured and estimated response was found to 
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be good for all the motion variables . One such typical 
matching for noise level of 5% is shown in Fig. 6. The 
results for all the four variables u, ^ , q and 9 show that the 
estimated response approaches the noisy measured response 
such that the random fluctuations of flight data due to 
presence of noise are quite evenly distributed on either 
sides of the estimated response. 

4. 3. 2.1 Results .for Case 2 

As mentioned earlier, the aircraft parameters from 
Table 9 were used as a priori fixed values while store 
parameters were estimated from the noisy response due to 
stores release. Following this procedure, store parameters 
were estimated for all four locations and for noise levels 
1%, 2% and 5%. A comparison of these results with that of 
Case 1 for corresponding store location ahd noise level is 
given in Tables 11 through 14. It is seen that the accuracy 
of parameter estimates in Case 2 has improved substantially, 
especially at higher noise levels (see results for- 5% noise 
level for all locations). However, at low noise levels, there 
is little to choose between results for Case 1 or Case 2, 
although the CR bounds for Case 2 are smaller as compared 
to those for Case 1. Further, in Case 2, fewer (2 or 3) 
iterations were required for convergence, making it computatio- 
nally economical. Matching of noisy measured response and 
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estimated response for one typical store location V2C2 
and noise level of 5% is shown in Pig* 7. The estimated 
response for all the motion variables looks to be matching 
well with the measured response except for fluctuations in the 
measured response due to noise* This observation is similar to 
that mentioned in $ 4.3.2. THe estimated response seems to 
a PP roac h the average smoothed response one would expect from 
the corresponding no— noise measured response. 

4.3.2. 2 Results for ca.se 3 

As in Case 2, here again aircraft parameters were kept 
fixed at values given in Table 9. Twenty measured responses 
with 20 different noise samples were analysed. Although 
results were obtained for all the four store locations, the 
results are presented for one typical location V1C2 and 
for noise levels varied from 1 % to 5%. Mean values and ct-^ 
for all the store parameters wet© obtained using Eqs. (3.1) 
and (3.2). The mean values and o~g of store parameters are 
shown in Table 15./This comparison is also shown in 
Table 15. A close look at the results for these two cases 
shows that the results from case 2 are marginally better 
than Case 3 for low noise levels (up to 2 %) whereas for 

/ It is of interest to compare the results for case 2 
and case 3. 
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higher noise levels , case 3 results were definitely superior 
to case 2, This can be explained easily based on the fact 
that the results for case 2 will be different depending on 
the chosen sample of noise to generate the measured responses 
It is, therefore, the results of case 3 which are represen- 
tative of all such possible results for case 2, in the sense 
of giving mean of all possible estimates from case 2. Thus, 
the comparison of Case 2 and Case 3 is of only academic 
interest. It may be emphasised again that the improvement 
obtained for Case 2 or Case 3 as compared to Case 1, is the 
most relevant fact to be observed because of its importance, 
to any practical application of the proposed scheme. Ratio 

Q ~ o 

lie between 0.74 and 1.2 indicating that estimated 

CR 

CR bounds are quite comparable to <7* as was mentioned 
earlier for case 2 also. 

4.3.2 .3 . Results for Case 4 

In this case, the aircraft parameters were fixed at 
values given in Table 10. The estimates for store parameters 
from noisy response due to store release are presented in 
Table 16 for noise levels of 1%, 5% and 10% for a typical 
store location of V2C1. The results for the above noise 
levels are also compared with those obtained for Case 2 in 
Table 16. It may be noticed that for all noise levels, 
the results for Case 4 are better than those for Case 2, 
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except for the parameter C whose accuracy deteriorates as 
noise level increases. However, even at 10% noise level, the 
matching between the estimated response and the measured 
response was good as shown in Fig. 8. As mentioned for Case 2 
(Fig. 7) , the estimated response again matches with what one 
would obtain from measured response after filtering out the 
noise. 

4 . 3 . 2 . 4 Results for Case 5 

In this case, aircraft parameters were fixed at their 
true values and only store parameters were estimated from 
noisy response due to stores release . Store parameters 
estimation was made for all the four store locations and 
noise was varied from 1% to 10%. A comparison among the 
results of case 2,4 and 5 is made for typical store location 
V2C1. Estimated values and CR bounds of all store parameters 
are given in Table 16. This table shows that the results 
obtained for Case 5 and Case 4 are almost identical, both ir. 
the estimated values as well as in their CR bounds. Thus, 
the observations made about the comparison between Case 4 
and Case 2 results also applies to comparison of results for 
Case 5 and Case 2. This case is of academic interest in the 
sense, that' it can not be used in practice since the true 
values of the aircraft parameters can never be postulated. 



However,, it may be mentioned that Case 5 represents the 
theoretical limit for Case 4 as far as expected accuracy of 
estimation is concerned. It also suggests the obvious that 
one should use the best possible estimates of aircraft 
parameter to analyse the noisy response due to stores release 
for obtaining the store parameters. 
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CHAPTER - V 


CONCLUSIONS 



7 \ method has been proposed for estimation of store 
parameters from measured response due to stores release. The 
Gauss-Newton method for parameter estimation was applied to 
a test case to show the applicability of the proposed scheme 
on simulated data. Results were obtained for varying noise 
levels in measured data, for different store locations and 
initial values of parameters. The effect of using Method 1 
and Method 2 on the accuracy of estimated parameters has also 
been presented. Based on the results presented/ we can draw 
The following important conclusions : 

(i) Application of the proposed scheme to no-noise response 

duo to stores release can estimate all parameters of 
aircraft as well as store in a single step following 
Method 1. Further/ accurate estimates with low CR bounas 
could also be obtained from the noisy response through 
method 1/ if the noise levels were low . 

(ii) For noisy response with high noise levels. Method 2 is 

recommended. The accuracy of store parameters can be 

increased by using the best available estimates of the 

aircraft parameters. To this purpose, many methods end 

10 

techniques have been reported in the literature where c 
one can improve the accuracy of the aircraft parameter 
estimates-these include design of proper flight maneuver 
choice of optimal input forms t instrumentation and 
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recording systems, and processing of data. Also, it must be 
mentioned that the proposed scheme for store parameter 
estimation can be used in conjunction with any one of the 
parameter estimation techniques (e.g. maximum likelihood, 
filter error methods, equation error methods, regression 
methods etc.) available in the literature. 

(iii) Por a given aircraft, the type of store and location 
will also influence the accuracy of store parameters 
estimates. Also, depending on the dynamic characteristics 
of the aircraft, a store release will dictate the resulting 
response of the aircraft following its release. Thus, any 
combination of aircraft and store that results in better 
dynamic response of the aircraft following stores release 
will yield relatively accurate estimates of store parameters 

(iv) A comparison of sample standard deviation (&~ ) of parameter 
estimates with Cramer-Rao bounds £ct^ r ) obtained in two 
different ways was carried out. Both cr s and cr CR were 
shown to be in reasonable agreement with each other. Thus, 
Cramer-Rao bounds are good practical measure for the 
sample standard deviation of parameter estimates and will 
provide a satisfactory indication of the reliability or 
confidence level of the estimated parameters. 
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***************** 

FILE: STR3 .FOR 
***************** 

0F° G SIBRE U FB0H T fLlGH° N TEST I, D»fs TU!,1N8i ' AER0DI ' ,SBIC COEFFICIENT? 

nm y P STUR| T PARftMPTrR^nt 1 V^i E FN£i E l§l £) At THEIR ESTIMATED VALUES 
REAL IYXA MA MS M |I ERS,CDS " CLS?CMS AKE ESTIMATED. this is method- * 

§*s 8 T* 6 ?| 6 ?^ ^ei^S?sg?f ^ ' 

i & & sf *< : lunvumv 



EXTERNAL UDOT , AL^DQT ,QDOT 
COMM0N/AI /XU , &TU, XAL£hA ,XU,X DELE 

epsffl ;s?;ss;:?K£^^x^g^feMEesf sfi ' ,,DEt ' E 

S 8 SS 8 B^^iifi;&§ ' 

QPEN(UNlT=I, D E VI CE='DSK',FlLE='RSP.DAT* f ) 

OPEN rUN!T=2f DEVICES *DSK';FILE=’TRSP. OUT') 

OPEN (UNIT=3, DEVICES # DSK%FILE='ERSP. OUT',) 

OPEN f UNlT=6 f DEVICE* 'DSK'lFlLE-fMRSP. OUT '5 

OPEN (UNlTs21,DEVlCE= f DSK f rFlLE='TPRSP. OUT') 

OPEN (UN%T*10 ,DEV |CE* f DSK* » FILE- 'RSP, OUT * ) 

I N = 1 •' ■ 

iTERsl 

READ (1,*)T0TTIM,TIMINC,U1,KC0UNT, NOISE, INMAX, ITRMAX 

ftPAH/l »nf.iu»rr n f .. ntlm n cr»tr» . 


READ f 1 , * ) ZU # Z ALPHA , ZALDOf # ZQ # ZDELE m 

read h , » ) mu » mttj,ma{pha,mtalpa ,MALDOT, MQ , MDELE 

READ(1»*5mS,CDS,CLS,CMS,SS,CS,D5,NS,XS 

IF(IN.NE.l)GO TO 8 


GO TO 10 
ISEED*ISEED+IN*57255 
GO TO 6 


— -I- 7 WRITE, 

WMTE (10, 11) 


INPUT VALUES 



E14 * 6/ 6X # * ZALDOTs * , E14.6/ 


^U»” ,C.l«*,D/OA# OA , AHurun- fE14.6/ 

V 6X; ' XQ* * , El 4 ; 6/6* , £XDELEsf ,Ei 4 _ 6 ) 

lao”li)fomM!flilNCfUl,KCOUNT,NOISEdNMAX#lTRMAX 
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WRITEU0.21 )MS,CDS,Cl.S,CMS,Ss f CS,DS,N5.XS 

r'?« A 4 C Hv' ,F r* i 4/5« '9Ri = ! 4^*4, 3X, ',CUS=',F7.4, JX. 'CMS=' , 

! ? S 2, 5 F5.2/3X, CS= »F5.2f 3X, f 0 S=' ,FSw5, 5x # 'NS= ' ,F4. 1 # 

1 jA| aps. #r / m 4 ) 

WRITE til) ,22 JTHEXAl.RH0,6 f SCF 

FnRHATu£,3HE|Al^F5.2;3X,'RHG=%F6.3,3X,'G=%F5.2, 
l A £ ? r O • 3 ) 

I^CtKCOUNI. E§.13. OR *(N0l$E. EQ.0})G0 TO 6 

‘ ~ ~ 71 4 J . ISEED 

, 5 / 

^10.6,5X,'ISEED= # # t i 2 > 

-initial ft.ight CONDITIONS — — " 

CONTINUE 

DATA THETA (1)/0,0/,UU 3/0. O/,QCl)/0.0/ 

DATA i^PHJCn/O.O/IPI/3.1416^ 

KTOT=98l’’"~"~~~~"~""™"~ 

NP=0,1/TIMINC 

AF=3.l416*DS*DS/4.0 

OB=0.5*RHO*U1*U1 

gs=s$F*yB 

DT=-V«S*CMS/MA) 

——estimation OF step inputs 

AX=(CDS*QS»AF*NS+MS*G*SIN(THETA1))/MA 

AZ= (C05*0S*SS*N5~MS*G*COS CTHETA1 3 3/MA „ „ # , T „v» 

AMx C (*CM 5 fCS«(XCC“XS)*CLS J *QS*SS^NStCDS*OS*AF*NS* CVHS"DT) ) /lYXA 

WRITEUO.183AX.AZf AM „ _ 

FORMAT! 5jc, »AX=* fF6.5f5Xf ?AZ=* fF8.5,5Xf •AM=*,F8.5) 

XDELE S AX .... 

ZDELE=A2i 
mdele=am 
GO TO 34 

REWIND? 21 ) 

READ (21 .*) (U2t CK) ,ALFA21 (K3 ,Q21(K3 ,TETA21 (K) ,K=1 f KTQT3 
irnowTssi ? ■ * *. * ■ 

m (KCOUNT.EQ.l ) ,0R. (KC00NT.EQ.2) ) JCOUNTsi 
Ei= 8 0S*AF*NS 
E2=QS*SS#NS 
DO 5 g 0 3 1 * JCQUNT 
T=Q.O. 



«• v m 

IFCCJ 

iFCil. _ 

if(j:eo,iS)Go 
xdelM.o 

ZDELE« 0.0 

8 B if 8 Sl * 0 

XDEL£« 1*0 
ZgELE*0.0 
MDEIjE*O #0 
GO TO 35 

xdele*ei/ma 
Z 5 E|iE» 0 j 0 (vH 5 -DT)/I y yA 


XDEltE^Oi© 

IPIrll 8 

XDEltE’OjiO 

S8H|2fa«.»)f«2/xm 
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GO TO 35 

xdele=o;o 

ZDEL£S(),0 
MDELE s -E2*CS/IYYA 
DO 50 K*1 , KTOT fi 
!£UKCOUNT,EqS}.0 

ifCuj.eq.i 3 . 5 ft im 
ifcui.sq.2i.or 
IFCtJ.EU. 35 .OR 
ifc CJ.eg.4).0R 
GO TO 806 
0 ELE =1 .0 

IFCKCOUWT.EQ.DGO TO 

if^kcount.e^dgo.to 



WP.EQ.DiGO TO 700 

• 0§.Cj.EQ.10nDgLE=ALFA21(K) 

D.FU.U))L>Ei»E=Q21£K) 
EU.123)DELE*1 .0 ’ 


( J i 


lF(MOl)(CK-w, 
KK=l+(K*l)/20 
PD(1 ,J,Kk5=UCK3/U1 

dl;WiBW K) 

sWi?n K,=THETStK) 


98 

too 


NP3 3110,108,110 


IFCMODC CK-1 3 .NP))1 10 -120,110 

alphas, $ 


F0RMATC4F12D) 
go To.no 



(K3 ,THETA£K) 


THETA CK) 


vr*. j »y viw -THETA(K) 

SOLUTION OF ‘EQUATIONS of motion by runge-kutta 

Sttfc 88SSH 

5A^A^2fV5i8g§IfR^itK^if^gS A<K> ' TBEXAt - KJ ' otK+1JJ 

T?T+fIMINC 


METHOD, 


3 , AlPHACK+1 ) ) 


CONTINUE 

F!SS!SSI:E8:i: 
xrt *£8SE5SftSSi8i 


IF ( ITER. 

ifcnqise 

N*49 
CALL ' “ 
CALL 
CALL 
CALL 
READ (2 
DO 60 
U(I)=U 


GO TO 1000 
GO TO 125 
GO TO 52 
I OF PSEUDO 
TO 54 


NORMAL NUMBER 


1 3 GO TO 54 


GGNOR i 
GGNORI 
GGNORi 
NOR I 

*5(1 

ggr' 

i: 


ISEED,N,Ri; 

iseed;n;r2; 


G 


.),Q(I3,THETAU),I=i,503 


'1 )*S|GC 1 ?*PERN 

>*5IGC23*PERN 
SRN 


O^I)*Qcl?+ii(I»l J*Iil?33*P! 

hthXbtmni fjaafcijrfi!j^ 8 Ri,.x.t, 5 P) 


.HU] 


! 0 i 8 ) 


WRIT|(6;6i 
FORMAT C4E2 
KC0UNT*3 
GO TO 33 
———•—FORMATION 
DO 13« 1*1,4 
DO 130 KK*£,50 
11*50* ( 1*1 )+KK 
DO 1|0 0X*i#3 

lF(J,GE,Dj*J+3 

IF(a,GE.55j*J+3 


OF MATRIX A 


mmWm * 
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SMU=SMU+(U«CI)-UE(I))**2 


IFCINMAX.NE.IUUJ to 132 
IF(ITER.$E,!tRMAX)GQ TO 132 
XFCNOISE.EQ.OJGO to 132 

0F CRAMER RAO BOUNDS — — 

ISGM A r 1 
IDP=5«> 

TFUsGMA.EQ.nGU TO 442 

PERU= 1 .0 

REWlNO(i) 

REvUND(6) 

WiP«f LPHAMt I JeQM(I) ,THETAM(I) ,1 = 1. 
READ(3,nCUECX},ADPHAE(I3,QE(I}>THETAEU3, 1=1,1 

SMt)=0.0 ’ 

DO 434 1 = 1 , IDP 

SMU=SMU+(UMCI)-UE(I))**2 

SIGU)=S0RT(SMU/FLOATCIDP)) 

$MA£»=0«0 
DO 436 1=1, IDP 

SMAL=SMAL+(ALPHAH( I ) -ALPHAE (I))**2 
SIGC2>=SQRTtSMAL/FL0AT(lDP)) 

smq=o;q 

DO 43B 1=1, 1DP 
SMQ=SMQ*CQMU)-QE(I))**2 ^ 
SlGC3)=SQRTCvSMQ/FD0AT(IDP)) 

SMTRsO.O 
DO 440 1=1 , IDP 

SMTH=SMTH+ { THETAM ( X ) -THETAECI ) ) **2 
SIGC4)=S0RT(SMTH/FL0ATCIDP)) 

Sl = (Ui*UD/C(PERN*SIG(l))**2> 

S2=1.0A (PERN*SIG(2) j**2) 

S3 = 1.0/UP£RN*SIG(3))**2) 
S4=l,0/t(PERN*SlGC4))**2) 

oo x t2®?(o I! ’ 200,I '* X *^ 

i-icHicKlSI'i^DlOGaNAL ELEMENTS of f matrix - 

DO 35 B 0=1,3 
SUMUsjO.O 
DO 350 1*1 ,50 
SUMO*SUHU + C A( 1 ,0 ) 3 **2 
Y1»SUMU*S1 

sumal»q;o 

DO 352 1*51,100, , 

SUMAE*SUMAE-KA(I,J))**2 

Y2=SUMAD*S2 

SUMQ*0.0 

DO 354 1=101,150 
SUH0aSUMO+( AU , J) )**2 
Y3»SUMQ*I3 
SUMT]H*0;0 

DEFfj)*mf.i+W+Y4 

D matr» 

DO 501 1*1,50 
DCX,I)=S1 • 

DO 502 1=51,100 
D C I , I ) =S2 
DO §03 1=101,150 


IDP) 

DP) 
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DCI,I)=S3 
DO S>04 1=151 , 200 
Oil# I )=S4 

lEiisKtoMSg U l 

no 262 i=t,IMX 
DO 262 d=l,LMX 

DO 264 0=1, UMX 
E ( 1 ,0 ) =ATC I »J)*D(U,J) 
COriTXHUE • 

——multiplication of 

DO 266 1=1 , LMX 
DO 266 K=l,LMX 
F(I,K)s0,0' 

DO 268 1=1 , LMX 
DO 268 K=1,LMX 
DO 268 0=1, JMX 

JVSIiffiM MHsfil 

CONTINUE 


263 

A MATRIX 


matrix 


MATRIX 


F0RMATC3X, 'DIAGONAL ELEMENTS OF F MATRIX* //CSX, 3E 18. 4)) 

— — — : inversion of f matrix — *—£*.. 

TUN Jt*3 r IFaIL=0 

CALL F01AAF(F,IA1,N,UNIT,IUNIT,WKSPCE,IFAIL) 

WRITE C5/DIFAIL 
DO 4|0 1=1,3 
CRBCi)=SORf CUNITC1.I)) 

WRITE 401 HCRBClSfl^lf 31 
IFCIsGMA.NE.DGO TO 444 ' 

FB0M MEASUREMENT N0ISEV/C5X, 

GO TO . 44f) J . . _ . _ 


1 /(5X.3E18.4) ) 


FROM ESTIMATED RESPONSE* 


* /■ V JA. J£,J 

ISGMA*ISGMA+ 

ifcisgma.gt. 

GO TQ 44§ 


REWINDC3) 

—FORMAT 


MATRIX B 
133 


M??5 E0 * 0,g ° 10 133 

iflilLD^ALPHAeiALPHAS 

iM 6 -OT 

8G"i8^» 

REWIND C 2 ) 1 50 

RKAP ALPHA? , Q2 r |H|TA2 
READ (3f *)U3# ALPHA3 , Q3 f THE I Ai 

lli$il!!« lAlP " A3 
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COMTI^UE 1 )=THETA2 " THETA3 

«R^wnr?? LUTIQN 0F LINEAR SIMULTANEOUS EQUATIONS 

M=2«o; IA=20O; IB-200? IDGTsb 
NA=3 
NB*1 

DO L I45 l ‘f = i ! R V A ' B ' m ' nA # n B,IA,IB, IDGT , WKAREA ? lER) 

xiia)=H(i!i ) 

;n;7rn-H Bl 2 A t^« OF STORE PARAMETERS — 

XOELE=CUS+XU , 1 ) 

CDS=Xl>ELE • 

ZDELE=CLS+XC2,1) 

CLS=ZOELE * 

MDELE=CMS+X13,1) 

cms=mdele ' 

«RITE(10.142)ITER 
FORMAT UOX, * ITERATION NO. M3) 

WRITER, *)£0S,CLs,CMS 
WRITE CIO, 146)CDS # CLS, CMS 

FORMAT(S| # JCDS= # ;F8.4,5X # 'CLSs%F8,4,9X^CMS=%r8.4) 

I X E.Ks* X ,i 1 

WRITE C5,*) ITER 
IfaTER.GT.lTRMAxJGO TO 200 
KCQUNT=2 
REWINDC21) 

REWlND(J) 

GO TO 17 

STORING OF ESTIMATED PARAMETERS — — — 

DD(1,1N)*CDS 
DD(2'lN)=CLS 
DD(3'lNI=CMS 
WRITE C5 j ♦) IN 
INsIN+1 

IF ( IN. GT. INMAX) GO TO 210 

rewindTi! 

REWIND(2 j 
GO TO 300 

CALCULATfo8*^ ,G MEAN ANf)°SAMPLE STANDARD DEVIATION 
DO 220 JN*1«3 
SUM(JN)*0,0 

88*?38ja5»t3B?48pc« f i giAX , 

DDM ( JN )*tsUM C JN 5 /FLOATC INMAX) 

SMTN(JN)=O f O 

STDVnC JN)»SQRT(SMTNt JN)/Fl«OAT(INMAX"l ) ) 

CuNTlNUK * ■ 

p0RMATf?it^MEA? M VALUfes N 0F'isTlMATED STORE PARAMETERS' 

IATION OF PARAMETER ESTIMATES*/ 

1 t7X # 3E20.4)) 

STOP 

Sitt signs 

CALL GGO0 
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CAUL MERFI 
END 

FUNCTION UnUT*" mm m m m mm>m m m m » m m. a .m m .m m m, m m m m 

FUNCTION UDOT (UX , ALPHAX /QX .ThFTAX ) 

COMMQ N/Al /XU , XTU J XALPhA XU *XDELE 

common/A2/mu,mtu'malpha1mt4i$a,haldot,mq m¥l E 

CCNMQN/A3/UU2U,§ALPHA»£Q,ZDFLE.7ALi)nT y ' MDELt 

q ’ 6UELE,ZAL00I 

CQMMON/A5/TIMINC 

UDOT=-G*THETAX+COS(THETAi) + CXU+XTU)^!JX+XALPHA*ALPHAX+XO»QX+XDEi.t, 

RETURN 

END 

* 9 $ 4m *W> mm mi *m «m m® «w> *** «* m® m* *m tmm ®m m* *m ®» HO mi on mat **** atm m& am mm %m mm me. fmm&imfm mi mm am tm mm « «w «n •» «» fn m 'm mb mb ,™. mb 

FUNCTION ALPDQT 

?8sssi9Si>ra:«i!M*;}jaiB!p“? 

S8358 taMi H0 ' MDELE 

COMMQ N/A4/G /DELE, THETA 1 
COMMQN/ A5 /TIM INC 

r p ^5^ESJSI)/?i5fi®f ,+zu * uxtzSLC ' HA * ALPHAX+UC!tul)f5X 

RETURN 

END 

HiNWONlooi 

«S£ I i8%taw;*«!ElfjIdl!H}?iQ.MB.E 

COMMON /A1 /XU, XTU, XALPHA /XU /XDELE 

COMMON/ A2/MU, MTU 'MALPHA.MTAlPA,MALDOT/MQ,MDELE 

C0MM0N/A3/Ut ,ZU/2ALPHA,Z0/ZDELE/ZALDOT 

COMMON/ A4/G/UELE, THETA 1 

COMM0N/A5/TIMINC 

ODD?! tMU+MTU)*UX+(MiLPfi^+ATALPAfiALPHAX+MALDOT*ALDOT+MQ*©X 

t ♦MDELE*DELE ’ 

return 

END_ 

subrqutine" , for" , Runge""kutta fourth order method" 

SUBROUTINE RUNKUtCFUH f VAR1 /VAR2 , VAR3 , VAR4/VAR5) 

COMMON/A5/H „ ' * 

Sl=sFUNCVARl,VAR2 / VAR3/VAR4 ] i 
S 2=FUN( VARl+0.5*H*Sl#VAR2/VARi/yAR4) 

S3=FUNC VARl+O.S^H^SS# VAR2f VAR3 , VAR4) 
S4sFUNCVARl+sl*H.VAR5*yAR3/VAR4) ' ' 

S«(Sl+S2*2,O+S34‘2,0TS4)/6iO 

VAR5«VAR1+S*H 

RETURN 

**♦ 


*SU****************************#**********************?***** 
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TABLE - 1 

Mass# moment of inertia, geometric characteristics and stability, 
control derivatives of test aircraft 

m a = 5000 kg ; I = 30400 kg_m 2 ; S = 16.42 m 2 ; 

yjr 'a 

b = 8.61 m ; C = 2.05 m ; U 1 = 137 * 5 m/sec 

© 1 = 0.0 rad g = 9.81 m/sec 2 ,* y> = 0.685 kg/m 3 

Mach = 0.44 ; = 4.2922 m 

Configuration clean ; Position of C.G. . : X_ _ = 25% of S 

C • G « 


Longitudinal stability and control derivatives 


C °u 

0.048 


0.0 

C °1 ■ • 

0.0375 

Cm l 

0.0 

C D 

0.41 

c 

-0.487 

C T 

X x 

u 

0.0 

c 

ra x 

-1.8 

c m 

\ 

0.0375 

C m 

q 

-4,4 

% 

0.0286 

c 

^e 

-0.67 


0.47 

G m T 

u 

0.0 

C L, 

Ps 

3*9 

C 

^rp 

‘‘'I 

0.0 

C T 

c< 

0.8 ' 

c 

m T 

0.0 

\ 

3.6 

C L,_ 

£>e 

0.32 
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TABLE - 2 


Arbitrary 

elevator input : 

time vs . deflection 


Time, sec 

$ e t rad 
(Positive 

Time, sec 

down) 

S e / rad 

0.0 

0.0 

1.1 

0.1 

0.1 

0.015 

1.2 

0.094 

0.2 

0.028 

1.3 

0.076 

0.3 

0.038 

1.4 

0.04 

0.4 

0.046 

1.5 

Q» 004 

0.5 

0.051 

1.6 

-0.024 

0.6 

0.057 

1.7 

-0.033 

0.7 

0.065 

1.8 

-0.026 

0.8 

0.086 

1.9 

-0.022 

0.9 

0.092 

2.0 

0.0 


TABLE - 3 


Standard deviation ( for various store locations 
corresponding to 5% noise level. 


CT ~ N 


Store Location 


V1C1 

V1C2 

V2C1 

V2C2 

G u • 

0.0195 

0.0115 

0.0200 

0.0140 

C £ 

0.0002 

0.0002 

0.0003 

0.0003 

O-q 

0.0003 

0.0003 

0.0003 

0.0003 

cr*» 

0.0007 

0.0005 

0.0008 

0.0006 



TABLE - 4 


True and 

Initial values 

of aircraft 

parameters for test case 


V 



SI. No. 

Aircraft 

Parameter 

True Value 

Initial Value 

1. 

-X 

u 

0.0074 

0.0037 

2. 

X * 

1*2759 

1.9140 ' 

3. 

6e 

0.0000 

0.0000 

4. 

- z u 

0.1498 

0.0749 

5. 

- z * 

83.7317 

41.8660 

6. 

“ Z q 

0.6975 

' 1.0463 

7. 

-z. 

6.8050 

3.4020 


h e 



8. 

M u 

0.0000 

0.0000 

9. 

- M * 

3.4918 

5.2377 

10. 

“ M q 

0.3314 

0.1657 

11- 

! 

© 

4.8040 

7.2060 



True and Initial values of store parameters and scale factor 
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-Rao bounds of parameters for location V2C2 : Method 



0274 .0.0294 0.0549 0.0588 0.1372 0.1470 



Mean values and c p* 0 for location V2C2 : Method 





8066 0*7517 0,8051 0,7617 0.8069 0.7531 
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CO P 
•H fO 
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TABLE - 9b 

Estimated values and CR bounds of aircraft parameters from a single 
measured response; input : elevator deflection. 
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